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In this paper we define the a posteriori probability W. This quan- 
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O ■ tity is introduced with the aim of suggesting a reliable interpretation 

of the actual threat posed by a newly discovered Near Earth Asteroid 
| (NEA), for which impacting orbital solutions (now commonly known 

as Virtual Impactors or Vis) can not be excluded on the basis of 
the available astrometric observations. The probability W is strictly 
related to the so-called background impact probability (that extrapo- 

■ lated through a statistical analysis from close encounters of the known 
. NEA population) and to the annual frequency with which the impact 
| monitoring systems (currently NEODyS-CLOMON at University of 

c/3 ' Pisa and SENTRY at NASA-JPL) find Vis among newly discovered 

asteroids. Of W we also provide a conservative estimate, which turns 
out to be of nearly the same order of the background impact prob- 
ability. This eventually says to us what we already know: the fact 
Q-r that nowadays monitoring systems frequently find Vis among newly 

discovered asteroids does not make NEAs more threatening than they 
have ever been. 



1 Introduction 

Soon after a new Near Earth Asteroid (NEA) is discovered, a preliminary 
orbit is computed using its positions in the sky over a suitable (minimal) in- 
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terval of time (astrometric observations). Like every physical measurement, 
astrometric ones are affected by errors which make the resulting orbit un- 
certain to some variable degree. Sophisticated mathematical and numerical 
tools are currently available to orbit computers which allow to propagate 
such measurement errors to the six orbital elements which identify the orbit 
of the asteroid. For this reason, the new NEA, soon after its discovery, is 
not represented by a single point in the 6-dimensional dynamical elements 
space; rather, it is represented by an uncertainty region, a 6-dimensional vol- 
ume with diffused contours. Obviously, the volume of this uncertainty region 
changes (usually shrinks) when new more observations become available and 
the orbit refines. 

Moreover, when the nominal orbit of the new NEA is geometrically close 
to the orbit of the Earth, and it shares some other peculiar orbital characteris- 
tics, it can happen that some orbital solutions which lead to a future collision 
of the asteroid with the Earth can not be excluded only on the basis of the 
available astrometric observations. Namely, orbital solutions which lead to a 
collision are inside the uncertainty region and they are fully compatible with 
the available astrometric observations and their errors. 

What is substantially done in these cases by the researchers, with various 
sophisticated techniques whose description is well beyond the scope of this 
paper (see Milani et ai, 2000; 2003), is to sample the uncertainty region 
according to a suitable frequency distribution (closely related to what is cur- 
rently known about error statistics) and then evaluate the relative probability 
that the "true" orbit of the asteroid is one of the collision ones. From now 
on we will refer to this probability with the symbol Vj. The collision orbits 
are nowadays commonly called Virtual Impactors (or Vis; for an exhaustive 
review see Milani et ai, 2003). 

Every time new more astrometric observations become available, the qual- 
ity of the asteroid orbit improves and the estimated impact probability V, 
is re-computed. Its value is almost always such that Vj <C 1 and during 
the phases of the orbital refinement it fluctuates, usually with a somewhat 
increasing trend 1 until it falls to zero, its most probable final value. 

Starting from 1999, some press announcements were made regarding as 
many newly discovered NEAs which were found to have non zero collision 

1 The reason of such increasing behavior is rather technical and it is essentially connected 
to the fact that uncertainty region shrinks with new more observations. 
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chances in the near future (given the highly chaotic dynamics involved in 
the multiple planetary close encounters, which are at the basis of the impact 
calculations, impact analysis procedures can safely cover only time spans 
of the order of a century). The computations were carried out mainly by 
two research groups, that at the University of Pisa and that at NASA-JPL. 
One of the first and, maybe, most famous of such cases was that of asteroid 
1999 AN 10 (for more information, see for example Milani et ai, 1999; 2003 
and http : //impact . arc .nasa. gov/news/1999/apr/21 .html; for a detailed 
historical account of these cases, see Chapman, 1999), being that of the 
asteroid 2002 NT7 the most recent one (up to this date). These objects 
obviously rose the somewhat alarmed attention of the public opinion and of 
the whole astronomical community for a while. Then, after the asteroids 
orbits were refined thanks to new more astrometric observations and the 
impact possibilities definitively ruled out, they have become again of purely 
academic interest. 

Currently, the only two existing automatized Vis monitoring systems, 
CLOMON 2 at University of Pisa and SENTRY 3 at NASA-JPL, find tens of 
newly discovered NEAs with Vis orbital solutions every year, and some with 
not so small impact probability (for a preliminary statistics of Vis detections 
see Tab. 1, more later). 

Given such past experience of public (and professional) reactions and 
given the current rate of Vis orbital solutions discoveries, some questions 
rise to the author's mind: how are Vis impact probabilities actually related 
to the real impact threat? How much threat should we reliably read in a Vis 
detection announcement? Equivalently, soon after the discovery of Vis or- 
bital solutions of a new asteroid, what is the probability that V approaches 
and eventually reaches the unity (within this paper we will use the com- 
pact notation "V — ► 1"), after the right amount of new more astrometric 
observations has become available? 

In this paper we give a statistical, a posteriori reading of Vis impact 
probabilities (which actually are in their very nature "deterministic", or, 
more properly, a priori) in order to provide an answer to such questions. 

2 http : //newt on . dm.unipi . it/neodys 
3 http : //neo . jpl . nasa . gov/risk/ 
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2 Statistical reading of Vi 



Soon after the discovery of Vis orbital solutions of a new asteroid, what is 
the probability of V« — » 1, after the right amount of new more astrometric 
observations has become available? It would seem quite obvious that such 
probability is simply Vi, according to the definition of Vj. But we believe that 
this is not the case. This is essentially because Vi, as we have said before, 
fluctuates every time new more astrometric observations become available 
and the computations are redone: which particular value should we consider 
for our needs? The value obtained with the second batch of astrometric 
observations following the discovery observations? Or the third? Or just 
the value of Vi calculated with the discovery observations? And, in latter 
case, if the discovery was made in another period of the year, would the 
calculated value of Vi have been the same? Actually, what we are asking 
is: only knowing that a newly discovered NEA exhibits some Vis orbital 
solutions, what is the probability that Vi will be equal to 1 at the end of 
the whole orbital refinement process? We believe that only a statistical, a 
posteriori analysis can give an acceptable and verifiable answer. 

Now, let us make the following thought experiment, only functional to 
the presentation of our point. Suppose that we are able to discover all as- 
teroids with absolute magnitude less than or equal to H which pass close to 
the Earth. Moreover, we reasonably suppose that every discovered impact- 
ing asteroid will show some Vis, with low Vi soon after the discovery and 
fluctuating with an increasing trend as soon as subsequent astrometric ob- 
servations become available. In other words, we are putting ourselves in the 
somewhat idealized situation where every impacting asteroid brighter than 
H is surely discovered and for it Vis monitoring systems surely spot some 
Vis soon after its discovery. 

Thus, we define the a posteriori probability of Vi — >■ 1, which could be also 
interpreted as a kind of "weight" of the Vis impact probability calculation 
(more on this later), as: 



W(<H)- n ^V 



P*(< H) m 
v(<H) T f Vi (<H)' U 

where n(< H) is the number of impacts of asteroids with absolute magnitude 
less than or equal to H and v(< H) is the number of asteroids with same 
absolute magnitude found among all the newly discovered NEAs to exhibit 
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Vis orbital solutions, both counted in the period of T years. Note that, 
according to what we have said at the beginning of this section, the number 
n{< H) is counted in the number v{< H), since we have assumed that every 
impacting asteroid is identified soon after its discovery as having some VI 
orbital solutions. 

Let us explain better the meaning of eq. (1). Within the hypotheses 
introduced above on the almost perfect NEAs discovery efficiency and Vis 
monitoring capabilities, we imagine to be able to wait for a very long period of 
years, T, and count the number n(< H) of asteroid impacts and the number 
v(< H) of Vis orbital solutions detected among all the discovered NEAs 
below the characteristic absolute magnitude H, within that period of time. 
Hence the fraction of these two numbers gives the a posteriori probability 
that a newly discovered asteroid, for which the Vis monitoring systems have 
spotted some Vis orbital solutions, is just that which will fall on the Earth. 
In the third member of eq. (1) we rewrite W in terms of the background 
annual impact frequency pi, namely that extrapolated through a statistical 
analysis from close encounters of the known NEA population (Morrison et 
al, 2003), and the annual frequency fy { of finding Vis among the newly 
discovered NEAs. A sketch of the average time between impacts (1/pi) is 
given in Fig. 1 as a function of the impactor's diameter. 

Before we have said that W can be interpreted as a kind of "weight" of 
the Vis impact probability calculation. This is because the greater is /v ; , the 
lesser is the a posteriori probability of V« — > 1, no matter how is the initial 
numerical value of V«. Namely, the higher is the frequency with which we 
find Vis among newly discovered asteroids (with respect to the background 
frequency of impacts, pi), the lesser is their a posteriori weight in expressing 
the threat of those particular newly discovered asteroids. This mechanism 
shares inevitable analogies with the well-known "crying wolf" experience, 
bad faith apart. 

Moreover, we can see that W is not directly related to the specific numer- 
ical value of Vj. Rather, it depends upon /y., namely the annual rate of Vis 
discovery, which, in turn, depends upon some observational characteristics. 
These are the annual number of NEA discoveries, the amount of astrometric 
observations available at discovery, the magnitude of astrometric errors and 
the conventions in their statistical treatment, as well as the observational 
geometry and orbital characteristics of the newly discovered asteroid. But, 
we guess, it is not easy to give an exact estimate of its value at the moment. 
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Average time between 
Impacts with Earth (years) 



Figure 1: An approximation of the average time between asteroidal impacts 
with the Earth (the reciprocal of the background annual impact frequency 
Pi) as a function of the impactor's diameter. We choose to use diame- 
ter, D, rather than absolute magnitude, H, since diameter is a more di- 
rect physical quantity. The mathematical relation between D and H is: 
log 10 D ~ log 10 1329 — y — \ log 10 £v, where pv is the albedo of the asteroid 
and it is usually assumed to be equal to 0.15 (Chesley et ai, 2002). 

A larger sample of Vis detections is necessary in order to better estimate fy. 
and thus the probability W. 

Anyway, having an idea of the total number of Vis detections found 
at every size between calendar years 2000 and 2001 (a sufficient but not 
complete archive of Vis detections could be found in the Observing Cam- 
paigns web pages of the Spaceguard Central Node at http: //spaceguard. 
rm. iasf . cnr . it), a conservative estimate of W cannot be too much differ- 
ent from that we get with /y. of the order of the unity. From Tab. 1 we see 
that the numerical value of /y ; varies between ~ 1 and ~ 10, in the reported 
range of absolute magnitudes H. For the sake of simplicity, in this paper we 
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choose to adopt /y 4 ~ 1 for all H, being confident of committing a justifi- 
able approximation, surely comparable with the uncertainty with which p« is 
currently known (at least within some ranges of absolute magnitude). With 
this reasonable assumption, we simply have W ~ p; L . 

Note that relaxing the optimistic assumptions on the almost perfect NEAs 
discovery efficiency and Vis monitoring capabilities makes fy v as approxi- 
mated with the aid of Tab. 1, even an underestimate. And consequently, it 
makes W ~ pi an overestimate. 

Therefore, the probability of V« — * 1 is nearly of the order of the back- 
ground impact probability, no matter how is V«'s specific, initial (fluctuating) 
numerical value. 

This result should not be a surprise since it simply re-states what we al- 
ready know: the actual impact threat of a newly discovered NEA exhibiting 
Vis orbital solutions is always the same estimated through the close encoun- 
ters statistics of known population. The fact that in the last few years many 
Vis orbital solutions have been detected among newly discovered NEAs ob- 
viously does not make NEAs more threatening than they have ever been. 

As a matter of fact, the annual rate of Vis detections, if compared with the 
background impact probability, suggests an order of magnitude of the weight 
Vis detections have in expressing the real threat of a newly discovered NEA 
with Vi ^ 0. 

3 Conclusions 

From the definition and the discussion of the a posteriori probability W done 
in this paper it follows that, rigorously speaking, the Vis impact probability 
Vi does not give the real expression of the actual impact threat posed by a 
newly discovered NEA exhibiting Vis orbital solutions. This is properly done 
by W, which is strictly related to the so-called background impact probabil- 
ity pi (that extrapolated through a statistical analysis from close encounters 
of the known NEA population) and to the annual frequency with which the 
impact monitoring systems (currently NEODyS-CLOMON at University of 
Pisa and SENTRY at NASA-JPL) find Vis among newly discovered aster- 
oids. Of W we also provide a conservative estimate, which turns out to be 
of nearly the same order of pi. 

All this might seem a bit paradoxical, given the definition of Vj. Yet, a 
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Table 1: Annual frequency of Vis detections below absolute magnitude H, 
estimated using the Spaceguard Central Node archive of Vis observing cam- 
paigns organized in the calendar years 2000 and 2001. During that period of 
time there were no Vis detections below H = 16. In the reported range of 
H, the numerical value of fy t varies between ~ 1 and ~ 10. 
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closer look at the definition of W shows that our conclusions are straightfor- 
ward and even obvious. 
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